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Astrophysical fluids may acquire non-zero electrical charge because of strong irradiation or charge 
separation in a magnetic field. In this case, electromagnetic and gravitational forces may act together 
and produce new equilibrium configurations, which are different from the uncharged ones. Following 
our previous studies of charged test particles and uncharged perfect fluid tori encircling compact ob- 
jects, wc introduce here a simple test model of a charged perfect fluid torus in strong gravitational 
and electromagnetic fields. In contrast to ideal iiia.giietoliydrodynarnic models, we consider here 
the opposite limit of negligible conductivity, where the charges are tied completely to the moving 
matter. This is an extreme limiting case which can provide a useful reference against which to com- 
pare subsequent more complicated astrophysically-motivated calculations. To clearly demonstrate 
the features of our model, we construct three-dimensional axisymmetric charged toroidal configura- 
tions around Reissner-Nordstrom black holes and compare them with equivalent configurations of 
electrically neutral tori. 



I. INTRODUCTION 

Equilibrium toroidal configurations of perfect fluid 
play an important role in studies of geometrically thick 
accretion discs around compact objects [I]. The isobaric 
surfaces also have toroidal topology and in order for ac- 
cretion to occur there must be a critical, marginally- 
closed isobaric surface with a cusp through which matter 
can outflow from the disc onto the compact object. In 
the following, we focus on black-hole systems and ignore 
self-gravity of the disc material. Shapes and properties 
of the tori, such as pressure and density profiles, arc then 
determined by the black-hole spacctimc geometry, an ap- 
propriately chosen rotation law (giving the distribution 
of specific angular momentum), and the fluid parameters. 

Perfect fluid tori in Schwarzschild and Kerr back- 
grounds were extensively discussed in the original funda- 
mental papers establishing this line of work [2, 3]. Later 
on, many studies appeared generalizing these models 
and including further details [4-10], describing tori also 
in Schwarzschild-de Sitter, Kerr-de Sitter and Reissner- 
Nordstrom-dc Sitter spacetimes, where presence of the 
so-called static radius [II] implies also the existence of 
tori with an outer cusp. 
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The material in accretion discs contains charged par- 
ticles (which may or may not be quasi-neutral in bulk) 
and the central black hole might also be charged. The 
charged or quasi-neutral fluid creates its own electromag- 
netic field which would then couple with that of the black 
hole, leading to a different and much more complicated 
description of the motion. Using the equations for the 
dynamics of the fluid and specifying its 'internal' proper- 
ties (conductivity, viscosity, equation of state, etc), one 
can solve the system so as to obtain profiles for the four- 
velocity, pressure, matter density and charge density [12]. 
However, the system of equations is rather complex, and 
in general requires the use of sophisticated numerical ap- 
proaches and codes, even if simplifying assumptions are 
made such as taking inflnite electrical conductivity (the 
limit of 'ideal magnetohydrodynamics'), no self-gravity, 
etc. On the other hand, some characteristic features of 
the motion of quasi-neutral or charged fluid, can also be 
studied relatively simply in a semi-analytic way [13 20]. 

The approximation of ideal magnetohydrodynamics 
(MHD) is reasonable in many astrophysically-rclevant 
situations involving fluids in motion [21]. However, there 
are other physical circumstances in which it is impor- 
tant to include the effects of flnite conductivity [22-24], 
and there the behavior becomes more complex, especially 
when strong gravitational and external magnetic flelds 
are also present. In order to address some of these effects 
in their mutual interplay, it can be useful to look also 
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at the opposite limit to that of ideal MHD: the Umit of 
neghgibly small conductivity. 

Here, we examine the problem of the interaction be- 
tween charged moving matter and the gravitational and 
electrostatic fields of the black hole, concentrating on an 
idealized situation which allows us to illustrate some oth- 
erwise very complicated effects. We consider a simple test 
model in which the matter is taken to be slightly charged 
and electrically non-conductive (dielectric), with the aim 
of providing an extreme reference case against which to 
compare subsequent more detailed calculations. We pro- 
ceed by first specifying a prescribed form for the angular 
momentum distribution, and then solving the dynamical 
equations to find the profiles of pressure, mass density 
and charge density. This approach can be seen as gen- 
eralizing the studies of uncharged perfect-fluid tori men- 
tioned above by adding the charge. Note that, through- 
out, our tori arc assumed to be composed of test fluids in 
which both the self-gravitational and sclf-clcctromagnetic 
fields are neglected. This gives a useful simplification, 
which is acceptable for weakly-charged, low-mass tori 
that have hardly any effect on the spacetime geometry or 
the ambient electromagnetic field. This approach helps 
us in building a semi-analytic model. In principle, in- 
clusion of self-gravity of the tori (following, e.g. [25-27]) 
and self-electromagnetic effects (see, e.g. [16]) would be 
possible, but that would enormously complicate the sit- 
uation. 

The charged dielectric perfcct-fliiid tori can also be 
seen as generalizing studies of charged test particles or- 
biting around black holes [28-36]. Various aspects of 
the charged test particle motion were also discussed in 
[37, 38], concerning the possibility of coUimated ejection 
along the axis of a rotating magnetized black hole, while 
investigation of stable off-equatorial lobes of charged par- 
ticles was discussed in [39-41]. In general, the motion of 
test particles is bounded within effective potential wells 
and this represents a model for a very dilute toroidal 
structure consisting of non-interacting particles; here, 
we add the non-electrical interaction between them, the 
pressure. Commonly, pressure leads to geometrically 
thick structures extending further out of the equatorial 
plane. 

In section II, we present the basic equations. In sec- 
tion III, we apply them to the case of a torus around 
a charged, non-rotating black hole described by the 
Reissner-Nordstrom metric. Nevertheless, the presented 
approach is suitable for a description of charged tori near 
to any models of compact objects with well-defined ge- 
ometry and electromagnetic field. We chose the Reissner- 
Nordstrom black hole because of its extremely clear ex- 
ternal electric field and geometry, given in an analytic 
form. For illustrative purposes we set the charge of the 
central black hole considerably exceeding astrophysically 
realistic values; our paper presents a toy model exhibiting 
the physical mechanism of mutual interaction between 
charged fluid and a black hole. In section IV, we dis- 
cuss the form of the isobaric surfaces for a torus with 



constant specific angular momentum composed of an un- 
charged barotropic perfect fluid. This is then extended 
to charged tori in section V, where wc also present a 
comparison between equivalent charged and uncharged 
cases. This work involves making a number of simplify- 
ing assumptions, and we discuss the nature and impact 
of these (including the zero-conductivity limit) in section 
VI. Section VII is the conclusion. Throughout the paper, 
we use the geometrical system of units (c = G = 1) and 
metric signature -t-2. 

II. BASIC EQUATIONS 

In general, the motion of charged perfect fluid is de- 
scribed by two sets of general relativistic MHD equations. 
These are the conservation laws and Maxwell's equations 

WpT°'^ = 0, (1) 
TJpF"^ = An J", (2) 

where the 4-current density J", which satisfies the con- 
tinuity equation 

V„ J« = 0, (3) 

can be expressed in terms of the charge density q, con- 
ductivity a and 4- velocity f/" of the fluid by using Ohm's 
law 

r =qU°' +aF°'I^Up, (4) 

with the electromagnetic tensor i^"^ being given in terms 
of the vector potential by F^p = Vq^^ — V pA^. This 
electromagnetic tensor describes the vacuum external 
electromagnetic field of the compact object (which per- 
vades the fluid), and also the internal electromagnetic 
fleld of the fluid itself, i.e., 

F'^P=F^i^ + F^^^. (5) 

The stress-energy tensor T'^^ consists of matter and elec- 
tromagnetic parts 

rpaP _ rpaP I rpafi /p^\ 

~ -'mat "T -'EM' V^J 

where 

Tmat = {e + pW'^U^+pg'^f', (7) 

^Et = ^ [f-^F^'' - \F,sF-^'g'^P^ . (8) 

Besides the pressure p and energy density e, the other 
fluid variables are the rest-mass density p and the spe- 
ciflc internal energy e = e/p — 1. The thermodynamical 
description of the fluid is specified by supplying an appro- 
priate equation of state p = p{e,q), which also involves 
the charge density of the fluid, describing the contribu- 
tion of the Coulomb interaction between the fluid parti- 
cles to the total pressure. 



3 



We build our model by considering a non-conductive 
{a = 0) charged test fluid (taken to be a perfect fluid) 
in an axially symmetric spacetime and use spherical 
polar coordinates {t,r,9,(j)). The fluid rotates in the 
0-direction with 4-velocity [/" = (f/*, f/'^, 0, 0), specific 
angular momentum £ = —U^/Ut and angular velocity 
(related to distant observers) fl = W^/U*, related by the 
formulae 



n 



igt t + gt4, 



t) = 



a 14 - Qua 00 

f-gu + 2lgt<p + c 



(9) 
(10) 



By writing out the covariant derivative of the electro- 
magnetic part of the stress-energy tensor (8) appearing 
the left hand side of the conservation law (1), moving it 
to the right hand side, and using the Maxwell equations 
(2), with 



a/3 _ 



INT 



47rJ° 



(11) 



we obtain the equation V/3T^^rp = F"^ Jp, with being 
the 4-current density due to the motion of the charged 
fluid torus. We are not here including the effects of the 
electromagnetic field generated by this 4-current: our tori 
are considered as being composed of 'test matter' from 
the elc(;troniagnetic point of view as well as from the 
gravitational one, i.e., Fj 
pap 



INT ^ -^EXT ^'^'^ '^^^ write 



-^EXT- Then wc get the master equation 

V7 T^"/3 _ 7" 
^^-'MAT ~ -^EXT-^/S- 



(12) 



Note that in this approach we do not need to solve 
Maxwell's equations (2), since the electromagnetic field 
is prescribed. Also, because of the non-conductivity, we 
have 



perfect fluid (sec, e.g., papers [2, 3]). In this case, when 
the angular momentum distribution I = i{r, 0) is chosen, 
a solution of the Euler equation for any barotropic fluid 
(having p = p{e)) can be derived from Boyer's condition 
[2, 3] 



dp 
P + e 



(16) 



W,^-W = ln(t/Oi,-ln(C/,)+ / j^^, (17) 

where the subscript 'in' refers to the inner edge of the 
torus in the equatorial plane. This condition enables us 
to straightforwardly determine the isobaric surfaces in 
the torus in terms of the equipotential surfaces of the 
'gravito-centrifugal' potential W{r,9): for equilibrium 
toroidal configurations composed of barotropic perfect 
fluid, the equipotential surfaces of W correspond to sur- 
faces of constant pressure (or energy density) in the fluid. 

When q ^ 0, the situation is more complicated and 
equations (14) are no longer integrablc for arbitrary 
q{r,0) = const and arbitrarily chosen I = i{r,9). It 
is then necessary either to specify I = £{r, 0) (with 
even £ = const being possible) and find an appropriate 
q = q{r, 6) which is consistent with that or, vice versa, 
to specify q = q{r, 6) (with even q = const being possi- 
ble) and find an appropriate £ = £{r, 9) . However, this 
is strictly 'necessary' only if the equation of state is pre- 
scribed. Otherwise, one coidd also absorb the constraint 
into that. The charged tori must clearly have distribu- 
tions of charge and angular momentum which satisfy the 
integrability condition. 



III. ISOBARIC SURFACES IN 
REISSNER-NORDSTROM GEOMETRY 



(13) 



The equations of motion (12) give two non-linear par- 
tial differential equations for the pressure p, whose pro- 
files we are wanting to find: 



drP = -ie+p)(dr ln(C/t)- 



deP = -{e + p)[de ln([/,) 



ndr£ 

1-Q£ 

nde£ 

i-n£ 



-qFraU" = n, 



qFga,U'' = T, 



(14) 



where TZ = TZ{r, 9) and T = T(r, 9). These equations are 
not integrable unless the integrability condition 



dell = drT 



(15) 



is satisfied, and this is therefore a requirement. 

The existence of a solution is guaranteed \i q = 0, 

so that the last terms in equations (14) vanish and we 
get the Euler equation describing a rotating uncharged 



The Reissner-Nordstrom metric, representing the 
space-time outside a charged, non-rotating black hole, 
provides a suitable mathematically-simple test example 
for experimenting with ideas before moving on to more 
complicated examples having direct astrophysical rele- 
vance. In the dimensionless (M = 1) Schwarzschild co- 
ordinates, the only free parameter in the line element of 
the Reissner-Nordstrom geometry 



r J 
-hr2(d6'^ -hsin^ edci)^) 



dr^ 

(18) 



is the dimensionless charge Q, which takes values \Q\ < 
1 for black- hole spacetimes and |Q| > 1 for naked- 
singularity spacctimcs. The locations of event hori- 
zons, corresponding to the pseudo-singularities of the 
geometry, are given by solutions of the equation A = 
— 2r + = 0. The ambient electric field is static 
and spherically symmetric, like the space-time, and is 
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described by the vector potential Aa = (^t,0, 0, 0), with 
the non-zero component being given by 



(19) 



In this background, the pressure equations (14) reduce 
to the form 



drP = -{P + e)(dr\n\Ut\- ] + q U'drAt 



deP = -{p + e) { de Inptl - 



1 - m 



where 



rsm 



0y/A 



(20) 



(21) 



To proceed further, it is then necessary to choose an 
equation of state. For an uncharged perfect fluid, a suit- 
able choice is to have a polytropic relation between the 
pressure and the rest-mass density 



p = Kp 



(22) 



with K and F being a polytropic constant and index. 
This widely-used relation is a convenient simple form 
which embodies conservation of entropy (as appropriate 
for a perfect fluid). For our calculations, we have used 
K = 10^^ and r = 2. This value of T is mathematically 
convenient for making analytic integrations and, while it 
is rather high for physical applications, the convenience 
makes its use consistent within the spirit of the present 
simple model. We have chosen a high value of k, be- 
cause electrostatic corrections to the equation of state 
then become negligible, so that we can continue to use 
this polytropic equation of state consistently even in the 
charged case (we comment further on this in section VI 
C). Moreover, we use values of p which are suSiciently low 
so that the medium is non-relativistic and the contribu- 
tion of the internal energy to the total energy density 
is then negligible as well, i.e., e ~ p. This approxima- 
tion is consistent also with the assumption of negligible 
self-electromagnetic-ficld. 

In order to find a solution for the pressure p, it is useful 
to rewrite equations (20) as equivalent equations for the 
density 



drP 



dep 



{np^ + p) (dr ln\Ut\- ^) - q WdrAt 
{Kp^ + p) (a,ln|[/,|-i^) 



-Tkp 



,r-i 



(23) 



which can be solved more easily than the ones for the 

pressure. On the other hand, in the uncharged case, the 
pressure profiles are determined from relation (16) which. 



for a polytropic equation of state, gives the following rel- 
atively simple analytic formula 



(24) 



(Note that this formula is valid only in the region where 

Win > W). 

In the next sections, we consider the commonly-used 
condition of constant specific angular momentum, £ = 
const [2-8, 8-10]. Tori with £ = const are particularly 
simple mathematically and are generally representative 
of those with a more general angular momentum profile, 
although some care needs to be taken when considering 
perturbations (these configurations are only marginally 
stable with respect to convective instability [42]). 



IV. UNCHARGED TORI, £ = const 

For a barotropic fluid, the isobaric surfaces coincide 
with the cquipotcntial surfaces of the potential W{r,6), 
which are given by the formula 



W{r,e) = In \Ut\= In 



r\/Asin9 



\J sin^ 



e-£'^A 



(25) 



One of the reality conditions here, A > 0, restricts the 
existence of equipotential surfaces to the stationary re- 
gion of the spacetime. The other one, r"* sin^ 9 — £^A > 0, 
evaluated in the equatorial plane {9 = tt/2), corresponds 
to the restriction that non-massless particles must move 
more slowly than photons 



< 



p2 



(r;g') 



r 

A' 



(26) 



where the function £'^^^{r;Q^) plays the role of the effec- 
tive potential governing photon geodesic motion in the 
equatorial plane (see [43] where the more general case 
with a non-zero cosmological constant is discussed). For 
the purpose of classification, we only need to consider 
positive values of ^ph('';<9^), i-e., we define £ph{f",Q'^) = 
The function £phi'r',Q'^) has one local cxtrcmum 
(a minimum) outside of the outer black-hole horizon, 

^ph,c(<9^), located at rph,c = ^(3 -|- \/9 - 8Q^), corre- 
sponding to the circular photon orbit in the equatorial 
plane. 

The character of the equipotential surfaces is well rep- 
resented by the behavior of the potential W{r, 9) in 
the equatorial plane, i.e., by the function W^/2(?') = 
W{r,9 = 7r/2). Since the orbits with vanishing gradient 
of W, i.e. those satisfying the conditions drW{r, 9) = 
deW{r, 9) = 0, correspond to loci with zero pressure gra- 
dients, the fluid has to follow geodesic motion there. The 
local cxtrcma of W are located only in the equatorial 
plane. Evaluating the necessary condition 9rW^/2(^) = 
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FIG. 1: Behavior of the function £c{r;Q'^) governing ex- 
trema of the potential W(r,9) in the equatorial plane of the 
Reissner- Nordstrom spacetirnc with Q = 0.1. For a fixed 
value of the angular momentum £, we can determine the po- 
sition of the potential maxima (smaller radius) and minima 
(larger radius). The potential W{r,9) is not defined in the 
shaded region limited by the function £ph{r;Q'^). In black- 
hole spacctimcs {Q'^ < 1), there is only the class of the mu- 
tual behavior of the functions £c{r; Q^) and ^ph(r; Q^), shown 
in this figure. The horizontal dashed lines denote the values 
of ^ms(Q^) = 3.670, ^„,b(Q^) = 3.995 and lph,c((?^) = 5.187 
for Q = 0.1, the value of Q being considered here, and the 
vertical solid line shows the position of the outer black-hole 
horizon. From the discussion of the behavior of W{r, 9) (see 
Fig. 2), it follows that equilibrium toroidal configurations of 
a barotropic fluid with £ = const exist only for £ > £m8(<5^)- 

0, these extrema are given by the condition 

i'=el{r;Q')^ ''^'^^^'\ (27) 

Again, for the purpose of classification, we only need to 
consider positive values of £c{r;Q'^). This function has 
one local minimum fms(Q^) corresponding to the spe- 
cific angular momentum of the marginally stable orbit. 
In addition to the limiting values ^ph,c(Q^) and £ms(Q^), 
there is also another one, inihiQ^)^ corresponding to the 
specific angular momentum of a particle moving along 
the marginally bound equatorial circular geodesic, deter- 
mined by the conditions drWj^/2{r) = and Wj^/2{r) = 0. 

Independently of the charge parameter \Q\ < 1 (we 
do not consider naked singularity spacetimes here) , there 
is only one type of mutual behavior for £ph(r;(5^) and 
^c(?';Q^)- We show this in Fig. 1, drawn for Q = 0.1 
which we take as a standard value in the following. (Note 
that wliilc this value for Q is rather high from an as- 
trophysical point of view [44], it still gives only very 
small deviations of the space-time away from that of the 
Schwarzschild metric; taking a value this high is useful for 
clearly demonstrating the effects which we are wanting 
to investigate.) 

Equilibrium toroidal configurations of a barotropic 
fluid with constant specific angular momentum exist 



only for £ > €ms(Q^), which takes values ranging from 
^ms(O) = 3.674, corresponding to the Schwarschild limit, 
to /ms(l) = 3.079, corresponding to the extreme Reissner- 
Nordstrom black hole limit. Moreover, an important fea- 
ture of constant specific angular momentum tori is that 
those which can form a cusp are limited by I < ImhiQ"^), 
which takes values ranging from Zmb(O) = 4 to /mb(l) = 
3.330. The limit ?mb(0) = 4 for the Schwarzschild case 
was also found for self-gravitating tori [26]. The radii 
given by the condition i = £c{r; Q"^) then correspond 
to motion of fluid elements along the unstable circular 
geodesic (smaller radius) and the stable one (larger ra- 
dius). The stable circular geodesic represents the 'center' 
of the torus (the potential W{r,9) has a local minimum 
there while the pressure is maximal there). The unstable 
circular geodesic marks a critical point (cusp), where the 
potential W{r,9) has a local maximum; the correspond- 
ing equipotential surface is self-crossing and is referred to 
as the 'critical surface'. As well as this critical equipoten- 
tial surface, there is also the characteristic null equipo- 
tential surface W{r, 9) = 0, which crosses the equatorial 
plane at infinity. 

The behavior of the potential W{r, 6) can be summa- 
rized in the following way: 

For £ G (0,£ms)) there are no extrema of the po- 
tential W, and there are no closed equipotential surfaces 
and no critical equipotential surface. The null equipo- 
tential surface is open towards the black hole (Fig. 2A). 

For £ = fms, there is one inflexion point of the po- 
tential W in the equatorial plane at which the critical 
surface has its critical point, corresponding to a ring. 
The null equipotential surface is open towards the black 
hole (Fig. 2B). 

For i e (^ms,^mb), there is a negative local maxi- 
mum and a negative local minimum of the potential W 
in the equatorial plane. In this case, closed equipotential 
surfaces exist which are bounded by the critical surface 
that self-crosses at the inner cusp. The null equipotential 
surfaces is open towards the black hole (Fig. 2C). 

For I = £mb! there is a zero local maximum and a 
negative local minimum of the potential W in the 
equatorial plane. The closed equipotential surfaces are 
bounded by the critical surface which coincides with the 
null equipotential surface (Fig. 2D). 

For £ e (^mb,^ph,c), there is a positive local maxi- 
mum and a negative local minimum of the potential 
W in the equatorial plane. The closed equipotential 
surfaces are bounded by the outer null equipotential 
surface. The critical surface is now open outwards 
away from the black hole, and self-crosses between the 
radii where the null surfaces cross the equatorial plane 
(Fig. 2E). 
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FIG. 2: Typical behavior of the potential W{r, 9) shown in terms of poloidal sections through the equipotential surfaces and 
the equatorial profile W^,r/2('") in a Rcissncr-Nordstrom spacctimc with parameter Q = 0.1. Taking progressively increasing 
values of the specific angular momentum £, samples are shown of 4+2 qualitatively different types of behavior of the potential, 
differing in the properties of the critical (dashed) equipotential surface and the null (thick) equipotential surface. 



For £ = iph.ci the potential W diverges at rph,c 
and the local maximum no longer exists. The negative 
local minimum of the potential W is still present. The 
closed equipotential surfaces are bounded by the outer 
null equipotential surface. The critical surface is no 
longer present. 

For £ > £ph.c, the only extremum of the potential 
W is the negative minimum. The closed equipotential 
surfaces are bounded by the outer null equipotential 
surface. There is no longer any critical surface, but there 
is a forbidden region for fluid elements with prescribed 
specific angular momentum, delimited by the radii 
satisfying the relation £ = tphir; Q^) (Fig. 2F). 



The behavior of the potential W{r, 9) is qualita- 
tively the same as in the pure Schwarzschild case [3]. 
The additional charge of the black hole Q here influences 
the values of £ms{Q^), ^mb(Q^) and ^ph,c(Q^)- 

Profiles of the pressure and mass-density can now be 
determined from relations (24) and (22). For doing this, 
it is necessary to choose relevant values for the parame- 
ters of the polytropic equation of state (22), for the spe- 
cific angular momentum and for the location of the inner 
edge of the torus. Here we present two examples, for tori 
with r = 2, «; = 10^2 ^nd i = 3.8 (Fig. 3): 

The first is the marginally bounded torus (the most 
extended closed torus), with its inner edge being at the 
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FIG. 3: Profiles of the potential W{r,9), pressure p(r,9) and rest-mass density p{r,6) of uncharged tori, shown in terms 
of poloidal sections through the equipotential, isobaric and iso-density surfaces, respectively, and their equatorial behavior 
W^7r/2(7'), Pir/2('') and p^/2(7'), for a Reissner-Nordstrom spacetime with Q = 0.1. Two examples are shown with different values 
for the radius at the inner edge of the torus: rin — rcusp = 4.544 in the upper figures and rin = 6 in the lower ones. The centers 
(pressure maxima) of both tori are located at rccnt = 8.388. The shapes of the equipotential, isobaric and iso-density surfaces 
coincide, with the values of the different quantities on them being related by equation (24) and the equation of state (22). We 
show poloidal sections through three such surfaces (A, B and C). The thick curve (surface C) marks the zero-pressure (and 
zero-density) surface which bounds the torus. The shaded regions above the profiles of W^/2(r) indicate the physically relevant 
parts of the profiles, delimited by the inner and outer edges of the tori. The specific angular momentum £ = 3.8 for both tori, 
and the polytropic parameters are F = 2 and k = 10^^ in each case. 



cusp. Tin — Tcusp- This is the most interesting case, since 
it can be used as a model for the inner parts of a thick 
accretion disk from which matter can flow in a standard 
way onto the black hole. To obtain the position of the 
cusp, we solve equation (27), which yields two real roots 
above the event horizon; n = 4.544 and rn = 8.388. 
The first root corresponds to the position of the unstable 
circular geodesic, i.e. to the cusp, while the second one 
corresponds to the stable geodesic, i.e to the pressure 
maximum: the 'center' of the torus. 

Our second example is a torus with its inner edge lo- 
cated at Tin > Tcusp- Here we simply chose the position of 
the inner edge to be at r^^ = 6- Note that the location of 
the inner edge must be chosen so as to be in between the 
cusp (n) and the center of the torus (ni), the positions 
of both of which are determined from the potential W 
independently of the choice of r-m ■ 



V. CHARGED TORI, £ = const 

In order to obtain the pressure or density profiles from 
equations (20) or (23), it is necessary to determine the 



charge density distribution q{r, 9) in the torus. This must 
satisfy the integrability condition (15). By expressing 
q(r, 0) in the form 

q{r,e)^q^p{r,e)k[r,e), (28) 

where is a constant and A:(r, 9) is a correction function, 
and using F = 2, we can rewrite equation (23) in the form 

drP = -^^[{Kp+l)dr\n\Ut\-qokU'drA^, 

dep = -^(^{np +l)de\n \Ut\y (29) 

Note that, as we express by relation (28), it is feasible 
to take the charge density distribution as being directly 
proportional to the rest-mass density distribution and to 
a correction function which represents variations in the 
charge per unit mass qok{r,9) through the torus, as we 
discuss in section VI A. From the mathematical point of 
view, the correction function plays the role of an 'inte- 
gration factor', which must be chosen so that equations 
(29) are integrable. 
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Now, due to the integrability condition (15), the cor- 
rection function k{r, 6) has to satisfy the relation 

2 sin 9{i'^A - sin^ O)d0k + 3k f A cos = 0, (30) 

which can be solved to give 

7 /^^^-r^sin^eX'/^ 



k{r, 0) 



PA - 



(31) 



where 7(r) is a function representing a constant of in- 
tegration over 9 for a given value of k in the equatorial 
plane, k{r, 9 = 7r/2) = 7(r). For the purposes of this pa- 
per, it is convenient to choose j{r) = 1. The charge den- 
sity distribution function with the correction function in 
the form (31) ensures the integrability of equations (29), 
and thus the existence of a unique solution for p(r, 9) and 
p{r,9). 

Integrating the second of equations (29) over the lati- 
tude, we obtain the following expression for the rest-mass 
density: 

, , 2'/^KC(r^sm''9-eAY^^ -V^ , , 
P{r,Q) = ^ • (32) 

The unknown function C(r), which depends only on the 
radial coordinate r, stands as a constant of this integra- 
tion. Its value can be determined by substituting the 
density formula (32) into the first density equation (29) 
and assuming the charge-density distribution according 
to relations (28) and (31). This leads to the ordinary 
differential equation 



IrAdrC (2r2 - 3r Q'^)C = 



AQqo 



2l/4K(r4 - £2A)3/4 



(33) 



Unfortunately, there is no analytic solution for equation 
(33), and so C(r) must be determined numerically. Since 
the torus is delimited by the zero-pressure (and zero- 
density) surface, we can find the necessary initial con- 
dition from the fact that p(rin) = 0. Prom the relation 
(32) we obtain 



C(rin) 



1 



2i/4«[rf„-£2(r2^-2rin + Q2)]i/4 



(34) 



As mentioned earlier, the most interesting configura- 
tion is the one delimited by the self-crossing zero isobaric 
(and iso-density) surface, i.e., the marginally bounded 
torus, which has its inner edge in the equatorial plane 
coincident with the cusp (rin = ''cusp)- We will be con- 
centrating on this type of torus from now on. The values 
of the specific angular momentum £ and the charges Q 
and qo then completely determine the shape of the torus 
and the positions of its center and of its inner and outer 
edges. At the cusp, the pressure and density vanish and 
have a saddle point (a minimum in the r-direction and a 



maximum in the 0-direction), in contrast with the cen- 
ter of the torus, where there is a local maximum. (In 
general, the density at the inner edge of the torus must 
be zero and when the inner edge is also a cusp, the den- 
sity profile has extrema there.) The location of the inner 
edge Tin, which needs to be known in order to evaluate 
the condition (34), can then be obtained from the first of 
equations (29) by setting drp\^_^ — 0, p(rin) = 0, and 
9 = 7r/2. This gives the following implicit expression for 



qoQrlyjA{ri-e-^A) + {Q^ - r,^)rt + fA^ = 0. (35) 

Note that when qo = 0, we get the relation {Q^ —rin)rf^ + 
i'^A'^ = 0, in agreement with the formula (27) derived for 
the case of the uncharged torus. 

To clearly illustrate how the charge on the torus affects 
its equilibrium structure, we constructed marginally- 
bounded tori with the same matter parameters and spe- 
cific angular momentum as in the uncharged case, i.e., 
we took K = 10^^ and i = 3.8, and we also considered 
the same charge of the central black hole Q = 0.1. Our 
sample tori are characterized by the parameters go = 0.4 
(positively charged) and qo = —0.4 (negatively charged), 
which, as we discuss in section VI A, correspond to the 
specific charges of the moving matter in the equatorial 
plane, due to our choice 7(r) = 1. 

For 50 = 0.4, equation (35) yields two real roots above 
the event horizon; n = 4.378 and rn = 9.104, where 
both of the roots correspond to circular orbits of charged 
test particles with the given specific charge 0.4. Choosing 
Tin = ^i, the numerical integration gives a regular thick 
charged torus. On the other hand, choosing Tin = rn, 
we get a degenerate torus (an infinitesimally thin ring) 
located just at r = ru- 

For qo — —0.4, equation (35) again yields two real roots 
above the event horizon; n = 4.767 and ru = 7.658, 
corresponding to circular orbits of charged test particles 
with the given specific charge —0.4. Again, the choice 
''in = Ti leads to the regular thick charged torus, while 
''in = '"11 gives the degenerate torus. 

In principle, one could also construct charged tori with 
ri < r-m < ru, as for the uncharged tori. However, this 
introduces further complications, and we focus here only 
on the more interesting critical (cusp) configurations. 

In Figs 4 and 6, we show the profiles of rest-mass den- 
sity and pressure. The related charge density distribu- 
tions q{r, 9) and correction functions k{r, 9) are shown in 
Figs 5 and 7. As can be seen from these, for the same Q 
and i, the positively charged tori are more extended than 
the uncharged ones while the negatively charged tori are 
less extended. Prom the right-hand panels of Pigs 4 and 
6, it can be seen that the densities in the more extended 
tori are larger than those in the less extended ones and 
so the masses of the more extended tori are clearly also 
larger. We come back to this in more detail in section VI 
B. 
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FIG. 4: Profiles of the pressure p{r, 9) and rest-mass density p(r, 6) for a positively charged torus (go = 0.4, I = 3.8) in a 
Reissner- Nordstrom spacetime with Q = 0.1, shown in terms of poloidal sections through the isobaric and iso-density surfaces 
(left), and in terms of their equatorial profiles p-n/2{i') and p^/2(f') (right). The torus terminates at r^^ — r^usp = 4.378 and 
Tout = 24.72, and its center is located at rccnt = 9.098. The dashed curves represent the zero-pressure surface for the equivalent 
uncharged case (left panel) and the profiles of density and pressure for the uncharged case (right panels). The isobaric and 
iso-density surfaces coincide (A, B and C), with the values of pressure and density on them being related by the equation of 
state. 



Q=0.1.f=3.8 Q=0.1,C=3.S. (7o=0-4 




5 10 15 20 25 5 10 15 20 25 

rxsin(^?) rxsm(^?) 



FIG. 5: Poloidal sections through the iso-contours for the correction function k{r, 9) (left) and the charge density q{r, 9) (right) 
for the same positively charged torus as in Fig. 4. 



VI. DISCUSSION 
A. Correction function and specific charge 

The dimensionless correction function fc(r, 9) has been 
introduced in equation (28) for convenience in calculating 



the torus configuration. It represents the variation with 
position of the specific electric charge (charge per unit 
mass), q — q()k{r,9). As shown in the previous section, 
setting fc = 1 so that the charge per unit mass is the same 
everywhere, does not give an equilibrium configuration: k 
must be allowed to vary so as to satisfy the integrabihty 
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Q=0.1./=3.8, qo^-OA 




rxsm(6) 



FIG. 6: Profiles of tlie pressure p{r, 9) and rest-mass density p(r, 6) for a negatively charged torus (go = —0.4, I — 3.8) in a 
Reissner-Nordstrom spacetime with Q = 0.1, shown in terms of poloidal sections through the isobaric and iso-density surfaces 
(left), and in terms of their equatorial profiles p,r/2ii') and Pti/2(.i') (right). The torus terminates at rin = ?"cusp = 4.767 and 
''out = 11.62, and its center is located at rcont = 7.660. The dashed curves represent the zero-pressure surface for the equivalent 
uncharged case (left panel) and the profiles of density and pressure for the uncharged case (right panels). The isobaric and 
iso-density surfaces coincide (A, B and C), with the values of pressure and density on them being related by the equation of 
state. 



Q=0.1,<'=3.8, 9o=-0.4 



A: k=\ 
B: A=0.999 
5 - C: A=0.996 
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A. g=-8>tl0" 
C: g=0 
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rxsin(^?) 



rxsin(f?) 



FIG. 7: Poloidal sections through the iso-contours for the correction function k[r, 9) (left) and the charge density g(r, 9) (right) 
for the same negatively charged torus as in Fig. 6. 



condition (15) and this requires the behavior (31). In 
the equatorial plane fc(r, 9 — 7r/2) — ■y{r), and the choice 
7(r) = 1 as a boundary condition is convenient for our 
present simplified model. From Figs (5) and (7), it can 
be seen that the required variations in fc(r, 9) away from 
this are actually very small (vifith the maximum being on 



the equatorial plane). 

Our choice of 7(r) — const (with the constant nor- 
mahsed to 1) corresponds to a case where the maximum 
of the charge density q{r,9) is located just at the center 
of the torus, where there is the maximum of the density 
p{r, 9), as can be seen from relation (28). Other choices of 
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7(r) could describe more physically relevant situations, 
but with the maximum of g(r, d) not necessarily being lo- 
cated at the center of the torus. For instance, by choosing 
7(r) = 1/r the specific net charge of the fluid in the torus 
grows in the equatorial plane from the outer edge to the 
inner edge, where it is maximal. Moreover, the maxi- 
mum of q{r, 9) is shifted from the center of the torus. Of 
course, at the inner edge, the net charge density g(r, 9) 
goes to zero together with the matter density p(r, 6). 

The values go = ±0.4 and 7 = 1, used for our rep- 
resentative cases, give a specific charge q in the torus 
around 10^** times smaller in magnitude than that for a 
proton, and so the medium can be thought of as having 
one particle in 10^^ with a net charge while the rest are 
neutral. We note that if wc decrease Q below our stan- 
dard value of 0.1 and simultaneously increase go in such a 
way that the product Qq remains unchanged, then we get 
essentially identical results. This is because the deviation 
of the space-time geometry away from Schwarzschild is 
extremely small for these values of Q, and so the rele- 
vant effect of the charge is almost entirely electromag- 
netic (depending on Qq) rather than having a significant 
gravitational contribution (depending just on Q). 



B. Total electric charge and mass of the torus 

In our model, we neglect the effects of the electromag- 
netic field generated by the torus, which is acceptable 
when this self-field is much weaker than the external elec- 
tromagnetic field associated with the central compact ob- 
ject. The total charge of the torus is given by 

Q = qy/^ dr dO d(t> 
Jv 

= Anqo I /' kp^f^dedr, (36) 

where 

9 or, = arcsm , (37) 
'° \^2K''r^C^ - I ) ^ ' 

is the function determining the upper boundary of the 
poloidal projection of the zero isodensity surface, g = 
—r^ siv? is the determinant of the metric tensor, and 
Tout is the radial position of the outer edge of the torus 
in the equatorial plane. For the representative tori which 
we are considering, with F = 2, k = 10^^, Q = 0.1 and 
(. = 3.8, we obtain the total charge on the torus as being 
Q = 2.49 X 10"-^-^ for go = 0.4 (positively charged) and 
Q = —4.50 X 10~^^ for go = —0.4 (negatively charged). 
Such small values, in comparison with the charge of the 
black hole, arc consistent with our neglect of the electro- 
magnetic field generated by the torus. 



The total rest-mass of the torus is given by 

M = py/^drd9d<l) 
Jv 

fTout f f 

= 4tt / p^/^dedr. (38) 

For the three representative tori which we are consid- 
ering, we obtain the total rest-mass as being M = 
6.24 X 10"" (positively charged), M = 1.13 x lO^^^ 
(negatively charged) and Ai = 9.64 x 10^^^ (imcharged). 
The ratio Q/Ai is then = ±0.4 in the positively and 
negatively charged cases, respectively, as clearly follows 
from the fact that the specific charge g — qok{r, 9) = 
±0.4fc(r, 9) = ±0.4, since k{r, 6) 1 throughout our tori. 

C. Parameters of the equation of state 

We set r = 2 for illustration purposes, since this choice 

simplifies the integration of the density equations (23). 
This value of F is inconveniently high for possible related 
astrophysical applications, but we stress that our model 
is an extremely siniplificxl one, purposely intended for 
investigating the behavior of a test case under extreme 
conditions (we are also taking a rather large value for the 
black- hole charge Q, only one sign of charge for particles 
in the torus and zero conductivity there). For such a test 
case, it would not be appropriate to introduce additional 
complications here in order to bring just one aspect of 
the model (the value of F) closer to astrophysical appli- 
cations. 

In general, electrostatic corrections should also be in- 
cluded in the equation of state, especially at higher mat- 
ter densities. However, this is not trivial to do (see, e.g. 
[45], for the electrostatic correction in the case of dusty 
plasmas). We have used a very high vahie of k, which 
enables us to neglect electrostatic corrections without in- 
consistency, because this enables us to demonstrate our 
approach more clearly. 

D. Distribution of specific angulcir momentum 

We have considered perfect-fluid tori with a prescribed 

specific angular momentum £ = —U^/Ut, which we set to 
be constant through the torus. In an uncharged case, L = 
and E = —Ut would be constants of motion connected 
with the assumed axial symmetry and stationarity of the 
spacetime. For a charged torus, the constants of motion 
are the generalized quantities L = + qA^ and E = 
-Ut - qAt. 

Note that the condition I = const is imposed for 

simplicity of the calculations; it is not essential for the 
method and can be relaxed. The Rayleigh criterion for 
linear stability against radial convection requires (. to be 
a non-decreasing function of the distance from the axis 
of rotation, and so £ = const uncharged tori are just on 
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the stability limit. For charged tori, the stability con- 
dition needs to be formulated in terms of a generalized 
quantity I = LjE [30, 46]. In the Reissner-Nordstrom 
electric field, the only non-zero component of the vector 
potential is At and one has 

l=^= ^ ^, , , (39) 

which reduces to £ for go = 0. Stability depends heavily 
on the specific charge distribution (i.e. on the behavior 
of the correction function fc(r, &)) and on the signs of the 
charges of the torus and the black hole. 



E. Approximation of negligible conductivity 

The assumption of high electrical conductivity of the 
medium is appropriate for many astrophysical plasmas 
with a high degree of ionization, and the ideal MHD 
framework can then be employed. Under the conditions 
of high conductivity and vanishing inertial effects of the 
plasma particles, the local electric field quickly becomes 
neutralized by rearranging the plasma flows (giving the 
conditions for the force- free approximation). A quasi- 
neutral medium then arises in which the volume density 
of net electric charge is negligible. 

However, there is an ongoing debate about the con- 
ditions that may lead to the presence of non-vanishing 
net charges, with an important role being played by elec- 
tric forces acting parallel to the magnetic field lines in 
the local co-moving frame. For example, a large-scale 
magnetic field may cause spatial separation of electric 
charges of different signs and their gradual accumulation 
in different parts of the system. Pulsar magnetospheres 
provide an example of such systems, with the charge sep- 
aration being caused by the dipole-type magnetic field of 
the neutron star [47, 48]. Black holes embedded in or- 
dered magnetic fields of external origin can also act in 
a similar way but then, for a low-density medium, the 
hydrodynamical description needs to be modified in or- 
der to describe the conditions of a coUisionless plasma 
(since the particle mean-free paths are then comparable 
with the characteristic length-scale of the system, given 
by the gravitational radius of the central black hole). 

One can imagine also another relevant scenario: a neu- 
tral fluid containing a few free charges, such as the case 
of dusty plasmas. In general, when charges feel an exter- 
nal electromagnetic field, they move generating a current, 
but if the fluid is dense enough and is highly collisional, 
the charges are less able to move and the conductivity 
becomes almost zero (see [49] for a discussion). Such a 
picture is actually compatible with our model since, as 
mentioned earlier, only a very small fraction of particles 
with net charge is required in order to give the parameter 
values used in our representative examples. 



F. Zero conductivity and consistency of the model 

The various limiting situations which we have been 
mentioning (hydrodynamical versus coUisionless plasma; 
infinite conductivity versus zero conductivity; self- 
gravitating matter versus test particles and fluids), are 
relevant under quite different circumstances and obvi- 
ously require different approaches. The approach which 
we have adopted in the present paper allows us to cap- 
ture the behavior of an idealized but non-trivial system 
where the fluid motion is governed by the combined ac- 
tion of a global (large-scale) electromagnetic field, the 
gravitational field of the central body, and pressure gra- 
dients operating within the fluid together with a non-zero 
electric charge distribution. 

The basic assumptions of our model are: 1) the fluid 
is a single-component test fluid (we ignore its self-gravity 
and its own electromagnetic field); 2) the fluid flow is 
stationary, with the 4-velocity having only time and az- 
imuthal components. If the conductivity a were non- 
zero, the second term on the right-hand side of the Ohm's 
law equation (4), which is proportional to tr, would give 
rise to a radial electric current unless there were a signif- 
icant self-fleld (contradicting the first basic assumption). 
Since our fluid is taken to be a single-species one, hav- 
ing a radial electric current would imply the existence 
of a radial mass current (contradicting the second basic 
assumption). Having cr = is therefore necessary for 
self-consistency. 

VII. CONCLUSIONS 

In this paper we have presented a model for a simple 

test case of an electrically-charged perfect-fluid torus ro- 
tating in strong gravitational and electromagnetic fields 
produced by a central compact object. Distributions 
of either specific angular momentum or charge density 
through the torus first need to be specified (we chose to 
specify the specific angular momentum distribution) and 
then pressure and density profiles can be calculated. An 
equation of state must be provided in order to close the 
set of equations. We have investigated the limiting case 
opposite to that of ideal magnetohydrodynamics, con- 
sidering a non- conductive (dielectric) perfect fluid rather 
than infinite conductivity as in ideal magnetohydrody- 
namics. Our case can describe a fluid in which bulk hy- 
drodynamic motion predominates over electromagnetic 
effects and the fluid has almost inflnite electric resistiv- 
ity. We are not here including the self-gravitational and 
self- electromagnetic fields of the matter in the torus, and 
so the treatment applies for low-mass, slightly-charged 
tori. 

For illustrating the model, we constructed both pos- 
itively and negatively charged barotropic tori, with a 
polytropic equation of state and constant specific angular 
momentum distribution, encircling a positively-charged 
Reissner-Nordstrom black hole. We compared the re- 
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suiting pressure and density profiles with an equivalent 
uncharged case, and also calculated the shapes of the 
tori. Taking the polytropic index F = 2, allows for the 
'pressure' equation to be integrated in a relatively simple 
way so as to give semi-analytic results. The large value 
for the polytropic constant k = 10^^ leads to tori where 
the Coulomb interaction between the charged particles of 
the fluid can be neglected in comparison with the stan- 
dard pressure due to the matter. The constructed tori are 
only slightly electrically charged in comparison with the 
charge of the black hole, thus generating a relatively weak 
electromagnetic field which can safely be neglected. It is 
striking that even with a very small value for the charge- 
to-mass ratio in the torus, significant effects are neverthe- 
less seen. However, it is necessary to stress out that the 
constructed tori carry the specific charge Q/M = ±0.4 
and orbit the black hole with Q/M = 0.1; since for any 
astrophysical body it would be practically impossible to 
maintain the specific charge Q/M > 10~^^ [44], the re- 
sults should be considered as illustrating samples only. 

The aim of this paper has been to introduce a new 
model of a dielectric charged torus encircling a charged 
compact object. We have proceeded by using a num- 
ber of simplifying assumptions: the compact object is a 
Reissner-Nordstrom black-hole; the torus is composed of 
test matter; the matter has a polytropic equation of state 
with prescribed F; the torus has constant specific angu- 



lar momentum, £ = const, and the specific charge is con- 
stant everywhere in the equatorial plane, 7 = 1. Despite 
the great simplifications coming from these assumptions, 
the scenario is still physically reasonable and non-trivial. 
Moreover, a large variety of degrees of freedom can be 
captured and the free parameters can be conveniently 
chosen in order to describe an astrophysically relevant 
situation. Of course, more complicated choices would 
then require more complicated calculations. Such calcu- 
lations are now in progress, but are beyond the scope of 
the present paper. 
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